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ABSTRACT 

Using the high resolution cosmological iV-body simulation MareNostrum Universe we 
study the orientation of shape and angular momentum of galaxy-size dark matter 
haloes around large voids. We find that haloes located on the shells of the largest 
cosmic voids have angular momenta that tend to be preferentially perpendicular to the 
direction that joins the centre of the halo and the centre of the void. This alignment has 
been found in spiral galaxies around voids using galaxy redshift surveys. We measure 
for the first time the strength of this alignment, showing how it falls off with increasing 
distance to the centre of the void. We also confirm the correlation between the intensity 
of this alignment and the halo mass. The analysis of the orientation of the halo main 
axes confirms the results of previous works. Moreover, we find a similar alignment for 
the baryonic matter inside dark matter haloes, which is much stronger in their inner 
parts. 

Key words: dark matter - galaxies: haloes - large-scale structure of Universe - 
cosmology: theory - methods: A^-body simulations - methods: statistical 



1 INTRODUCTION 

Understanding how the angular momentum is distributed 
over the collapsed structures in the Universe and its con- 
nection with the large-scale structure distribution of galax- 
ies at present is a key issue to determine its origin (see 
|Primacl3l2003l . for a review). Furthermore, angular momen- 
tum is fundamental in order to explain the structure and 
dynamics of galaxy discs and hence it is a basic ingredient 
to determine their evolution and morphology. The currently 
accepted paradigm for explaining the origin of the angu- 
lar m omentu m is t he Tidal Torque Theory l|Peebleslll969l : 
iDoroshkc vichi Il970l : [Whitc 1984), hereafter TTT. In this 
framework, the angular momentum is transferred to the pro- 
tohalo by the coupling of the quadrupolar inertia tensor and 
the tidal tensor during the linear regime of growth of density 
fluctuations. According to TTT, the epoch of turnaround 
marks a major change in the amount of angular momen- 
tum transferred, as this tidal torquing becomes inefficient 
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due to the small size of the halo after it collapses. How- 
ever the theory breaks at this epoch as non-linear effects 
are expected to turn up together with merger and accretion 
processes. In particular, this causes the dilution of the TTT 
prediction regarding the spin direction (|Porciani et al.|[2002l : 
[Pcirani ct al. 2004). 

None the less, recent observational studies have found 
a preferential orientation of galaxy discs with respect 
to their surrounding larg e-scale structures (|Navarro et al.l 
l2004l : iTruiillo et aLlT2006l ). Such an aUgnment is thus ex- 
pected to be a relic of the alignment already present be- 
fore the epoch of the turnaround. These works differ both 
in the metho dology and the sample used. For instance, 
iNavarro et al.l (2004) used catalogues of nearby galaxies and 
found an excess of edge-on spiral galaxies highly inclined rel- 
ative to the supergalactic plane. However this result has two 
major drawbacks: the sample is limited to galaxies in the 
vicinity of the Milky Way, which is not very representative 
of the overall large-scale structure; and the orientation of 
the large-scale sheetlike distributions like the supergalactic 
plane is blurred due to redshift distortion at high distances 
(whic h can only be determ ined by redshift). On the other 
hand, ITruiillo et al.l l|2006l ) used the largest galaxy redshift 
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surveys up to date (2dFGRS and SDSS) to get the sam- 
ple, and searched for large voids in these surveys to identify 
spatial orientations not relying on the detection of planes. 
They found that spiral galaxies located on the shells of the 
largest cosmic voids have rotation axes that lie preferentially 
on the surface of the void. As this paper has been written, 
it has also been found a correlation between the orientation 
of galaxy discs a nd local tidal sh ears using the Two Mass 
Redshift Survey (|Lee fc Erdogdu,2007 l. 

In addition, Af-body cosmological simulations have 
nowadays enough resolution to study whether dark matter 
haloes also show this alignment. This result is expected as it 
is commonly assumed that both dark and baryonic matter 
shared a similar evolution during early epochs and gained 
the same specifi c angular momentum b efore the formation 
of the disc (e.g. iFall fc Efstathioulll980l '). However, even in 
the case of a positive signal, it is expected to be very weak as 
dark matter haloes are more affected by torques from neigh- 
bouring haloes than baryonic matter concentrated in discs. 
This has recently lead to numerous studies with a variety 
of methods to select a homogeneous sample of dark matter 
haloes with the aim of searching a counterpart of the obser- 
vational results. For example, in lAragon-Calvo et al.l (|2007| ) 
a multiscale filter is developed in order to distinguish haloes 
belonging to filaments or walls. This code infers the mor- 
phology of the structure surrounding a halo by calculating 
the relations between the eigenvalues of the hessia n of the 
density field. In a dynamical approach, iHahn et al.l (|2007al ) 
use the number of positive eigenvalues of the hessian of the 
gravitational potential to classify the environment in which a 
halo resides. Both papers conclude that haloes in walls have 
spin vectors that tend to lie in the plane of the host wall, but 
haloes in filaments show only a weak trend for their angular 
momenta to be aligned with the filament direction. How- 
ever, the method proposed bv lTruiillo et al. I liooe) provides 
a very clean way to characterize the orientation of the large- 
scale distribution (i.e. the reference direction to measure any 
kind of alignment) by using the centre of a large void. This 
method, besides its robustness in observations with respect 
to redshift distortion, owes its convenience to the well known 
fact that the direction of maximum compression (that of the 
smallest eigenvalue of the deformation tensor) turns out to 
be very well correlated with the direction to the centre of 
the void. 

The study of this type of alignments is specially inter- 
esting in the context of gravitational lensing. A key assump- 
tion of this method is that the observed galaxy ellipticity 
correlations come only from the distortion caused by the 
gravitatio nal shear, with all the intrinsic terms being negli- 
gible fe.g. iHevmans et al.ll2006h . Thus, exploring these cos- 
mological alignments may prove to be useful in order to esti- 
mate the contamination from this source in the lensing sig- 
nal. The importance of angular momentum is also remark- 
able in semi-analytical modelling of galaxy formation (e.g. 
lAvila-Reese et al]|l998l ) . The properties of the galactic discs 
are usually expected to be related to those of the dark mat - 
ter haloes in which they are embedded (e.g. IMo et al.|[l998l ). 
However, regarding angul ar momentum this assumpti on has 
been recently questioned jP'Onghia fc Navarrci|l2007l ). 

In this p aper, mainly motivat ed by the observational re- 
sult found bv lTruiillo et al. we use a high-resolution 
cosmological simulation to study the alignment of shape 



and angular momentum of galaxy-size dark matter haloes 
in the shells of large voids. Our si mulation pro v ides en ough 
statistics to a ssess the results in iPatiri et al.l (|2006al ) and 
iBrunino et al.l (|2007,) and for the first time to measure the 
strength of the alignment signal of the angular momentum 
in these haloes with out performing a pre vious selection of 
them, as opposed to lBrunino et al.l l|2007l ). 

This work is organized as follows. In Section [2] we de- 
scribe the cosmological simulation and the methodology. 
The results are presented in Section O We summarize and 
discuss our conclusions in Section |4l 



2 MARENOSTRUM UNIVERSE 

SIMULATIONS AND METHODOLOGY 

The cosmological simulation we have used for this st udy is 
the MareNostrum Universe (|Gottlober fc Yepes^l2007^ . This 
non-radiative SPH simulation employs 1024^ dark matter 
particles of mass 8.24 x 10®/i~^Mq and 1024^ gas particles 
of mass 1.45 X lO^/i"^M0 m a box of 500ft" ^ Mpc on a 
side. The cosmological model is ACDM with cosmological 
parameters = 0.7, Q,m ~ 0.3, Q,bar = 0.045, h = 0.7, 
(jg = 0.9 and a slope n = 1 for the initial power spectrum. 
The evolution of the initial conditions was p erformed by us - 
ing the TREEPM+SPH code GADGET-2 (|Springe]||200"5l ). 
The spatial force resolution was set to an equivalent Plum- 
mer gravitational softening of 15 comoving kpc, and 
the SPH smoothing length was set to the 40*'' neighbour to 
each particle. The long range gravitational force calculation 
was done by the Particle-Mesh algorithm using EFT in an 
homogeneous mesh of 1024'^ elements. 

In this work we use both dark matter and gas particle 
distributions of this simulation. In order to test the efi'ects 
of the baryonic component in the results of our analysis, we 
have also used the results of an N-body, dark matter only 
version of the MareNostrum Universe. This simulation has 
exactly the initial conditions and parameters than the SPH 
simulation but with 1024^ dark matter particles only. 

Dark matter haloes are found using a parallel hierar- 
chical friends-of-friends algorithm based on the minimum 
spanning tree for the particle distribution with linking pa- 
rameter b = 0.17, which corresponds to overdensity A = 330 
with respect to mean ma tter density. Using equation 2 in 
iBailin fc Steinmet j l|2005h . and taking into account that the 
mean axis ratios for galaxy-size haloes are c/b ~ 0.79 and 
b/a ~ 0.75, we find that about 65 particles are enough to 
determine main axes orientations to within 15°. As this is 
a reasonable resolution, we are restricted to use haloes with 
M > 5.4 X lO"/i"^M0, so we will choose our lower mass 
limit above this value. There are 881,861 haloes in the en- 
tire simulation box above this mass at z = 0, which provides 
very good statistics for our purposes. 

The detection o f voids is made using the HB void finder 
(jPatiri et al.l l2006al '). This code searches for the maximal 
non-overlapping spheres which are larger than a given radius 
and empty of objects below a given mass. Those spheres are 
our voids. We are interested only in large voids, i.e. those 
with -Rvoid > 10h~^ Mpc, which characterize the orienta- 
tion of the large-scale distribution at the position of the 
halo. We searched for voids defined by haloes with masses 
larger than 7.95 x 10^^/i~^Mq, whose axes are resolved to 
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within 13° . The mass limit has been chosen to get the same 
number density ol haloes as the number de nsity of galax- 
ies fou nd in the observed galaxy samples of iTruiillo et al.l 
l|2006l ) (n = 5 X 10~^h^ Mpc~^, which corresponds to galax- 
ies brighter than Mtj —19.4 -I- 5 log h). With this defining 
mass we found a total of 3047 voids, with a median radius 
of 11.17/i~^ Mpc. For further comparison with the observa- 
tional result, we look for Milky Way-sized haloes in a shell of 
4h~^ Mpc around the surface of each void. We thus remove 
all the haloes with mass below the defining mass of the voids 
7.95 X 10^^ft~^MQ (which is greater than our resolution 
limit), and those with mass exceeding 7.95 x lO^^/i'^M©, 
so that they are still galaxy-size. Considering all these con- 
straints, we are finally left with 88,426 haloes in the shells 
of voids and in this range of mass for our analysis. 

For every halo we calculate its total angular momentum. 



L = ^mrfc X Vk 



(1) 



where rk has its origin in the centre of mass of the halo and 
the sum is over every dark matter particle identified by the 
halo finder as belonging to the halo. The orientation of the 
angular momentum with the centre of the void is given by 
the 9 angle derived from: 
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where R is the vector linking the centre of mass of the halo 
with the centre of the void, R = rhaio — f void. In the same 
way we use an analogous equation to calculate the angle 
between R and the major (a), middle (b) and minor (c) halo 
axis, i.e. the three main axes resulting of the diagonalization 
of the 'inertia' tensor of haloes in the shells of the voids (the 
eigenvalues will be referred to as a, b and c). This tensor is 
defined as: 



E 



(3) 



where the coordinates Xi,k represent the i-th component 
of the position vector of the fc-th particle, measured with 
respect to the centre o f mass of the halo. T his sum is again 
over halo particles (see lPatiri et al.ir2006bl '). 



3 RESULTS 

The main results are summarized in Fig. [T] and Table [1] We 
show in Fig. [T] the probability density distribution of the 
orientations of the three main axes and angular momentum 
of haloes in the shells of voids (using a spherical shell of 4/i~^ 
Mpc around the surface of the void), with respect to the 
direction to the centre of the void. The angles are calculated 
using equation [2] and the analogous ones for the main axes 
correspondingly. 

We also plot an analytical fit of the probability den- 
sity distribution (Betancort-Rijo & Trujillo, in preparation) 
given by the following equation: 

pdfj. 



P(^)d/i oc 



fi = cos 9 



(4) 



(H-(p2_ 1)^2)3/2' 

where p is a free parameter, related to the ratio between the 
dispersion in the radial component of the angular momen- 
tum vector, and its dispersion in the transverse component. 



Although this equation was proposed to describe how the 
orientation of angular momentum is distributed, it proves 
to be a very good fit for the distribution of the orientations 
of the halo main axes. 

The distribution of an isotropic distribution is com- 
puted and also shown in Fig. [T] If these vectors have no 
particular orientation (our null hypothesis), then the ex- 
pected probability distribution is a sine function. We use 
this fact to represent the deviation of the probability distri- 
butions with respect to the isotropic case, so that the null 
hypothesis corresponds to the zero level in this plot. 

In order to quantify the strength of these alignments, 
we use the paramet er p in Eq. |4j as previously done by 
iBrunino et al.l ()2007h . The null hypothesis corresponds to 
p = 1, so that is uniformly distributed. Lower and higher 
values of p correspond to a preferential orientation parallel 
and perpendicular to the centre of the void respectively. A 
good approximation for this parameter in the case p ~ 1 
comes from the relation: 



(cos 9) 



(5) 



where the average is taken over all the halo sample in the 
shells of voids. 

The results presented in Table [1] are for the three main 
axes and angular momentum of haloes in the shell of 4/i~^ 
Mpc around the surface of the void. The first column shows 
the result of averaging the cosine values of the angles con- 
sidered here for our sample of haloes. The second column 
shows the result of the statistical test performed here, i.e. 
the average of the cosine test. In the case of isotropic orienta- 
tions, the average of the cosine is normally distributed with 
mean 0.5 and variance 1/ (12A^), where A'^ is the number of 
haloes. Therefore, the bigger is the departure from the mean 
and the bigger the statistics, the higher the significance of 
rejecting the null hypothesis, = VT2iV((cosS) — 0.5). 
The third column represents the strength of the alignment 
using the aforementioned parametrization, estimated from 
a Levenberg-Marquardt fit to our data using Eq. |4l The 
fourth column also represents the value of p but this time 
it is estimated using Eq. [5] The error bars in this parame- 
ter correspond to the 68% confidence level. The last column 
shows that the value of the reduced is indeed near to 
unity, which proves the goodness of the fit. It is important 
to remark that the statistical test only quantifies the proba- 
bility of isotropy in these orientations; the magnitude of the 
alignment is instead measured by the parameter p. 

We find for this sample of haloes that the minor axis is 
preferentially aligned with the direction to the centre of the 
void, and the major axis is preferentially aligned perpendic- 
ular to it. Middle axis also shows an alignment orthogonal 
to the radial direction, but much weak er than previous one . 
This confir ms the r esults reported by Patiri et al.l l|2006lj ) 
and Brun ino et al.l l|2007h . The same conclusion holds from 
the analysis of the gas particles embedded in dark matter 
haloes, although in this case middle axis does not show any 
significant alignment. 

In order to characterize the complete statistical infor- 
mation about the spatial orientation of these haloes, we 
studied the distribution of the Euler angles. These angles 
link two orthogonal systems defined at location of each halo 
in the shell of a void: the one formed by the direction to the 
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Figure 1. Probability density distribution (measured with respect to the case of isotropic orientations) for the angles between the 
angular momentum of dark matter haloes in the shells of voids, the three principal axes of these haloes, and the direction to the centre 
of the void. The error bar on each bin represents the Poissonian error. The dashed line shows a fit of the data using Eq. |4] and the dotted 
line represents the null hypothesis. Taking haloes in the shells of spheres with randomly distributed centres makes the results agree with 
isotropy. 



Table 1. The strength of the alignment of the three main axes and angular momentum vector of haloes (and baryons embedded therein) 
in a shell of 4/i~^ Mpc around large voids, with respect to the direction to the centre of the void. First column shows the average of the 
cosine of the angles between these vectors and the radial direction R. Second and third columns show the value of the parameter p from 
a fit to Eq. |4]and using Eq. [5] respectively. The fourth column shows the goodness of this fit. The last column represents the result of 
the statistical test for the rejection of the null hypothesis in terms of A^o- (see text for details). 



axis 


(cos e) 


Pfit 


^estimated 


X^-iVd.o.t. 

V^'Vd.o.f. 




dark matter 


major 
minor 
middle 
L 


0.4739 
0.5306 
0.4945 
0.4926 


1.110 ±0.004 
0.884 ± 0.003 
1.023 ±0.004 
1.029 ±0.004 


1.110 ±0.004 
0.885 ± 0.003 
1.022 ± 0.004 
1.030 ± 0.004 


+0.06 
+0.63 
-0.04 
-0.69 


-26.9 
+31.6 
-5.7 
-7.6 


gas 


major 
minor 
middle 
L 


0.4639 
0.5288 
0.5042 
0.4922 


1.154 ±0.012 
0.891 ± 0.009 
0.980 ±0.010 
1.030 ±0.011 


1.156 ± 0.013 
0.891 ±0.010 
0.983 ± 0.011 
1.032 ± 0.011 


-0.56 
-0.36 
-1.15 
+3.20 


-13.1 
+10.4 
+1.5 
-2.8 



centre of the void and any two perpendicular directions that 
lie in the shell (X, "F, ii), and the one formed by the three 
main axes derived from the diagonalization of the inertia 
tensor of the halo (a,6, c). The second coordinate system 
can be obtained from the first by means of a set of three 



rotations: a rotation of angle about It, a rotation of angle 
about the rotated X axis (which places the Tl axis in the 
direction of a), and finally a rotation of angle t/; about this 
a axis, which places the other two in the direction of h and 
c. It is important to remark that lo is actually the angle 
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between the major axis and the direction to the centre of 
the void, so we aheady know its distribution. Moreover, due 
to the rotational symmetry about the direction R (there is 
no privileged point in the surface of a sphere), we also know 
that 4> is uniformly distributed. Thus, for a complete descrip- 
tion of the spatial orientation of our haloes we need to find 
out the distribution of the ip angle. Obviously if the shape 
of the haloes were not aligned with any particular direction, 
Ip would be uniformly distributed. Nevertheless we find that 
the distribution of this angle is instead well described by 
(Betancort-Rijo & Trujillo, in preparation) 

P (ip) dip ^ - {I + f3 cos {2'tP)) dtp (6) 

TV 

where /3 is a free parameter. A fit of our data gives (3 — 
— 0.093±0.005. However, if we select a sample of those haloes 
whose major axis forms an angle with R above a given value 
oJo, this parameter seems to rise. This is an indication that 
in fact, the distribution of ip is not independent of this angle, 
and the complete description of the shape alignment will be 
given by a more general distribution P{ij,tp). 

But our most remarkable result is that we find a pre- 
viously undetected alignment of angular momentum, whose 
orientation is preferentially perpendicular to the direction to 
the centre of the void (but not necessarily aligned with the 
major axis). This is confirmed by the average of the cosine 
and a Kolmogorov-Smirnov test. Although this alignment 
is small, both statistical tests agree in rejecting the null- 
hypothesis with a probability corresponding to a significance 
level of more than 7a. 

Previous studies have not found such an orientation of 
the angular momentum in the shells of voids, mainly be- 
cause of the poor halo statistics due to the small simulation 
box. For instance, the halo sample in IPatiri et al.1 (|2006bl l 
consisted of only 1729 objects with which they managed to 
find p — 0.96 ± 0.04. Despite this value is compatible with 
our result, it points that angular momentum tends to be 
aligned with the direction to the centre of the void, which is 
the opposite of what we find. B ox size may also have been 
the reason why iHevmans et al.l ((2006) could not reject the 
null hypothesis (causing the error bars to be about a fac- 
tor of two larger than those shown in our Fig. [1} , although 
other factors may have had some influence on their conclu- 
sions, like the small number of particles in their haloes or 
their es timation of error bar s via bootstrap method. The 
case of H runmo et al.1 (|2007h is specially surprising, since 
having the same box size, the resolution of Millennium sim- 
ulation largely exceeds that of this paper. Yet, they still find 
that the distribution of the angular momentum is compati- 
ble with a random orientation. Taking a convenient subsam- 
ple of haloes they found a significant alignment, but only 
in one out of the seven radial bins considered. This sub- 
sample is obtained selecting haloes with a disc-dominated 
ga laxy at their c entre, using semianalytic galaxy catalogue 
of ICroton et al.l ()2006l ). The galaxy must be brighter than 
Mk < -23 and have a bulge to total ratio < B/T < 0.4. 
On the contrary, it is noteworthy that no preselection was 
performed for the results presented here. 

The likely explanation for this disagreement is that their 
haloes were identified with the same halo finder but with 
a linking length which corresponds to an overdensity 900 
times the mean matter density, to focus on the core prop- 



erties of the haloes. Here, we compute main axes and an- 
gular momentum of haloes detected with a linking length 
corresponding to 330 times the mean matter density (see 
Section [2]). For comparison, we show in Table [2] the results 
of the same analysis as in Table [T] but using a catalogue 
of haloes in our simulation box defined by an overdensity 
eight times higher. Although the main axes show a stronger 
alignment in the inner parts of the halo, it is clear that the 
alignment of the angular momentum with the direction to 
the centre of the void is lost when we focus in the inner 
parts. The contribution to the angular momentum from the 
particles in the outer shells of a halo is then fundamental 
for this conclusion. It is quite remarkable that the gas in 
the inner part does not show this alignment either, despite 
the fact that both major and minor axes show an enhanced 
alignment compared to dark matter haloes. 

As a further test of our results we repeated the analysis 
to obtain the strength of the alignments, by taking haloes 
in the shells of fake voids (i.e. spheres with randomly dis- 
tributed centres throughout the simulation box). The results 
for both shape and angular momentum alignment are then 
compatible with the null hypothesis. To test the effect of 
baryons in the alignments of dark matter haloes we have 
analysed the dark matter only version of the MareNostrum 
Universe. The results we obtain for the strength of these 
alignments are compatible to those listed in the upper (dark 
matter) part of Tableland Table [2] 

We also studied the dependence of the strength of these 
alignments when moving outside the surface of the void. The 
results are shown in Tabled In this table we show the value 
of the parameter p estimated using Eq.Ualong with their er- 
ror bars, in several different shells surroundin g the void. The 
shells are selected for direct comparison with iBrunino et al.l 
l|2007l ). It is clear from the table that the strength of the 
alignments dilutes as we move outside the surface of the 
void. The inner shell show a weaker alignment of the angu- 
lar momentum compared to that work. Instead, this strength 
declines slower here as we advance to t he outer shells. This i s 
in contradiction with the statement in iBrunino et"aLl (|2007l ) 
suggesting that taking a wide shell can mask this alignment. 
Interestingly, we find that the alignment of the angular mo- 
mentum is still significant even at few times 7?void. 

The minimum void size i?void > 10/i"^ Mpc was chosen 
i n ord er to mimic the observational result of lTruiillo et al.l 
(|200d ). and it seems reasonable as it is of the order of the 
scale of nonlinearity at present, and hence haloes around 
voids with this size are embedded into quasi-linear struc- 
tures which remember the original shape. However, we are 
in an excellent position to study whether this value is good 
enough for our analysis. In Fig. [2]we show the dependence of 
the strength of the alignments on the choice of the void size, 
in the range of 5/i"^Mpc < _Rvoid < 15/i"^Mpc. The width 
of the bins is chosen to be Aiivoid = 2h~^ Mpc. All the 
haloes in a shell 0.4_Rvoid in width around the surface of the 
void are taken into account. We find that the alignment de- 
creases considerably with the radius of the void, fading away 
for -Rvoid < 5/i~^ Mpc. On the other hand, large voids are 
rare and the poor statistics of objects around them increases 
the uncertainty on the parameter p. Therefore, a strong and 
significant alignment can be found using voids of intermedi- 
ate size. It is worthy noting that if we use the same procedure 
as above, i.e. all the voids larger than a given minimum size 
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Table 2. Same as Table [T] but the haloes are defined here as enclosing eight times the virial overdensity. 



axis 


(cos e) 


Pfit 


Pcstimatcd 


X^-iVd.o.t. 
V2'Vd.o.f. 




dark matter 


major 
minor 
middle 
L 


0.4702 
0.5362 
0.4925 
0.5003 


1.125 ±0.005 
0.869 ± 0.004 
1.030 ±0.005 
1.000 ±0.005 


1.127 ±0.005 
0.865 ± 0.004 
1.030 ± 0.005 
0.999 ±0.005 


-0.23 
+1.00 
-1.42 
-0.70 


-25.2 
+30.7 
-6.3 
0.2 


gas 


major 
minor 
middle 
L 


0.4476 
0.5549 
0.4926 
0.5015 


1.232 ±0.020 
0.793 ±0.013 
1.029 ±0.017 
0.988 ±0.016 


1.234 ±0.021 
0.802 ±0.014 
1.030 ±0.018 
0.994 ±0.017 


-0.65 
+2.54 
+0.95 
-0.19 


-12.3 
+12.8 
-1.7 
+0.4 



Table 3. The strength of dark matter halo alignments measured from parameter p at Icr confidence level using different shells around the 
surface of the void. Recall that the null hypothesis (isotropic orientations) corresponds to p = 1. Lower and higher values of p correspond 
to a preferential orientation parallel and perpendicular to the centre of the void respectively. 



Shell 


Major axis 


Middle axis 


Minor axis 


Angular momentum 


Number of 


/?void Units 


P 


P 


P 


P 


haloes 


1.00 <R< 1.05 


1.138 ± 0.010 


1.022 ± 0.009 


0.865 ± 0.007 


1.046 ± 0.009 


16406 


1.05 < R < 1.10 


1.149 ± 0.013 


1.008 ± 0.011 


0.874 ±0.010 


1.024 ±0.012 


9437 


1.10 <R< 1.20 


1.115 ± 0.008 


1.019 ± 0.007 


0.883 ± 0.006 


1.021 ± 0.007 


24329 


1.20 <R< 1.40 


1.085 ± 0.005 


1.027 ± 0.004 


0.901 ± 0.004 


1.024 ±0.004 


66576 


1.40 <R< 1.80 


1.052 ± 0.003 


1.006 ± 0.003 


0.946 ± 0.002 


1.013 ±0.003 


187565 


1.80 <R< 2.60 


1.025 ± 0.002 


1.008 ± 0.002 


0.969 ± 0.002 


1.009 ± 0.002 


460592 


2.60 <R< 3.20 


1.013 ± 0.001 


1.001 ±0.001 


0.985 ± 0.001 


1.006 ± 0.001 


555621 



are included into the analysis, the highest signal-to-noise 
ratio is located around < Gh~^ Mpc. As shown in Fig. |31 
this ratio declines very steeply with the choice of minimum 
void size because the decrease in the statistics (and hence 
the increase in the uncertainty Ap) is much higher than the 
increase in the strength of the alignment. This result is of 
interest in observational studies, but it is important to re- 
mark that the position of this maximum could be shifted 
when using redshift space instead of real space. 

Some studies have suggested that voids are not gener- 
ally spherical, so that alignment detections assuming spher- 
ical voids should therefo re be seen as residual of genuine 
physical alignments (e.g. lAragon-Calvo et al]|2007l ). How- 
ever, using an ellipsoidal void catalogue for the MareNos- 
trum Universe we find similar results for the orientation of 
angular momentum and main axes for the haloes in the shells 
of those non-spherical voids (see Table |4}. Moreover, we find 
no significant difference between the direction to the centre 
of the void and the direction perpendicular to the ellipsoid, 
taking n = (R ■ a/a? ,R ■ b/b^ , R ■ c/c^) as the reference 
direction for angle measurements. Dark matter haloes se- 
lected here are located in a shell of 4/i~^ Mpc thick around 
the maximal sphere enclosed into the ellipsoid. In order to 
study more in depth the relevance of void shape, we studied 
the probability distribution function of the Euler angle (j). In 
the case of spherical voids, this angle is uniformly distributed 
due to the rotational symmetry in the plane tangential to 
the sphere at the position of the halo. On the other hand, 
this symmetry is broken in the case of ellipsoidal voids, so 



Table 4. Values for the p parameter for the alignment of the 
three main axes and angular momentum vector with the direc- 
tion to the centre of the void for haloes in a shell of Ah~^ Mpc 
using two different void catalogues in the same simulation. For 
the ellipsoidal void catalogue the results for both the radial and 
the normal direction are shown. 



axis 


spherical void catalogue 


ellipsoidal void catalogue 




p (using R= n) 


p (using Ft) 


p (using n) 


major 


1.110 ±0.004 


1.112 ±0.004 


1.112 ± 0.004 


minor 


0.884 ±0.003 


0.886 ± 0.004 


0.887 ± 0.003 


middle 


1.023 ±0.004 


1.021 ± 0.004 


1.021 ± 0.004 


L 


1.029 ±0.004 


1.028 ±0.004 


1.029 ± 0.004 



the distribution of the </> angle might be different. However 
taking into account that 

. V ■ i - {v ■ R){i ■ R) 



^l-{vRf^l~{i-R) 



(7) 



where i is defined in the direction of the major axis of the 
void and v in the direction of the major axis of the halo, 
we find {<j}) I {11/2) = 0.5001 ± 0.0009 in agreement with a 
uniform distribution. 

There also have been recent claims of dependence on 
the alignment of the angular momentum with the mass of 
the halo in filaments and walls of the large- scale structure 
(jAraeon-Calvo et al ]|2007l : lHahn" et al ]|2007al ). In agreement 
with these papers we find a slight trend of increasing align- 
ment strength with halo mass for our sample. Using a shell 
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Figure 2. Tlie strengtli of tlie alignment as a function of the void 
size, measured over bins of 2h~^ Mpc wide in void radius. The 
solid, long-dashed and short-dashed lines represent the major, 
middle and minor axes respectively. Dotted line represents the 
angular momentum, which has been slightly shifted to the right 
for the sake of clarity. Error bars show the la uncertainty in the 
value of the parameter p. The dash-dotted line represents the null 
hypothesis. 



Figure 3. The signal-to-noise ratio in the alignments of haloes 
around voids as a function of the minimum void size. This ratio is 
estimated by the deviation of the parameter p from the isotropic 
case p = 1, over the uncertainty in this parameter. As in Fig. |2] 
dotted line represents the angular momentum and the solid, long- 
dashed and short-dashed lines represent the major, middle and 
minor axes respectively. 



of 4/i ^ Mpc the haloes near the resolution limit imposed, 
with Af = 5 X 10"-8 x 10"/i"^Mq, have a strength value 
of p ~ 1.015. On the other hand, haloes with masses in the 
upper limit of the sample A/ = 8 x 10^^-2 x 10"/i"^Mq 
show evidence of a stronger alignment with p ~ 1.060. 



4 DISCUSSION AND CONCLUSION 

In this paper we have analysed the orientation of main axes 
and angular momentum of dark matter haloes in the shells 
of large voids using the MareNostrum Universe simulation. 
The spins of these haloes are preferentially perpendicular 
to the direction to the centre of the void. The possibil- 
ity of isotropic orientations has been rejected at 7a con- 
fidence level. This result is in qualitative agreement with 
and is likely to be related to the observational result of 
iTruiillo et al.l (|2006l ) for spiral galaxies using large galaxy 
surveys. The alignment of dark matter collapsed structures 
with the large-scale structure finds a natural explanation in 
the context of the Tidal Torque Theory. The different com- 
ponents of the angular momentum vector present a disper- 
sion which is not isotropic: there is a marked asymmetry be- 
tween the dispersion in the direction of maximum compres- 
sion as compared to that in the direction of maximum ex- 
pansion (Betancort-Rijo & Trujillo, in preparation). There- 
fore, very low density regions such as large voids are suit- 
able to study this effect as they are surrounded by a region 
with a large tidal field, causing the anisotropy in the shear 
tensor we are measuring in this work. It is worthy noting 



that this effect is not due to accretion of matter along fila- 
ments. Large voids arc surrounded by a complex filamentary 
network of matter (Colborg ct al. 2005) wh ich could repre- 
sent a preferential direction of infall (but see lVitvitska et al] 
l2002h . This accretion generates angular momentum pointing 
perpendicular to the filament in which the haloes are em- 
bedded, which could be either the direction to the centre of 
the void but could also be perpendicular to it, resulting in 
no preferential orientation of angular momentum. Therefore, 
the anisotropic dispersion in the angular momentum vector 
seems more plausible in order to explain of this effect. 

This is the first time that this alignment has been found 
in cosmological simulations using haloes in the shells of voids 
with no preselection of the haloes, in contrast to previous 
results. The poor halo statistics in the simulation box, the 
small number of particles in haloes, and even an overes- 
timation of the error bars are li kely to be the caus es of 
the non-detecti o n of t h is effect in Patiri et al.l (|2006lJ ) and 
iHevmans etall (120061 ). iBrunino et al.1 l|2007h did not suc- 



ceed in finding this alignment either, and concluded that 
the distribution of angular momentum vectors is compatible 
with isotropy for their total dark matter halo sample. This 
is somewhat surprising since they used the Millennium sim- 
ulation which exceeds the resolution of the MareNostrum 
Universe simulation used here, having both the same box 
size. However, we have checked that the uncertainties in the 
orientation of the vectors used here are small enough to en- 
sure the veracity of these results. We recall that the haloes 
used here are chosen to be well resolved and their directions 
are determined to within an angle of 13°. Indeed, it can be 
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proven that this small indetermination in the angle measure- 
ment makes the uncertainty in the parameter p be negligible. 
It is very likely that the linking length (and hence the over- 
density) used to identify the haloes in the simulation box 
makes the results to differ as we have shown in Section [S] 
The choice of an overdensity 900 times the mean matter den- 
sity (as opposed to our selection of overdensity 330) allowed 
them to focus only in the inner parts of the halo, but this 
may have prevented them from the detection of the align- 
ment of angular momentum we find. This suggests that the 
particles in the outer shells of the halo are determining in 
the measurement of the strength of this alignment. 

None the less, the claim of iBrunino et all (|2007D that 
they found a significant alignment at the surface of the void, 
is based only in the choice of a subsample based on a semian- 
alytic galaxy catalogue. As we have mentioned above, they 
do not find any preferential orientation in their total sample. 
Moreover, we find that this alignment is significant even at 
few times -Rvoid- The signal decreases with distance to the 
void centre, which is expected for several reasons: the disper- 
sion in the different components of the angular momentum 
vector becomes isotropic with increasing distance to the cen- 
tre of the void. This is due to the less anisotropic tidal field 
in regions with higher density. Besides, the local density in- 
creases at larger distance from the centre of the void, and 
hence the interaction with neighbouring haloes is stronger, 
diluting the alignment from the initial torquing. 

The choice of the minimum size of the void -Rvoid > 
10h~^ Mpc has proven to be reasonable because it is a good 
compromise between strength of the alignment and statis- 
tics. However, we have found that the choice i?void > 6/i^^ 
Mpc is even better for the detection of alignment of dark 
matter haloes. When including smaller voids, which are more 
abundant, the statistics of haloes in the sample is highly im- 
proved. On the other hand, the alignment of these haloes is 
much weaker compared to the strength of the haloes around 
the largest voids: haloes around smaller voids are separated 
by shorter distances, and thus they are more influenced by 
interactions. Thus, the best choice is non-trivial and requires 
from the large amount of data given by cosmological sim- 
ulations. Prospective observational studies should flnd this 
scale useful to get high signal-to-noise ratios in the detection 
of these alignments. 

Our main result that galaxy-size dark matter haloes 
around voids have spins th at lie in the shells of voi ds, 
is in agreement with that of lAragon-Calvo et al.1 (|2007l ) if 
haloes in the shells of void which form large- scale sheets 
are ma inly responsible for this effect. However, iHahn et al.l 
l|2007bl ) found a stronger alignment of haloes in sheets than 
largely exceed the one presented here. Nevertheless, haloes 
around voids are located in a less dense environment than 
actual sheet haloes, which may have influence on their ori- 
entation. 

Shape alignment of haloes in the surface o f large voids 
is furth er conflrmed, as shown previously by iPatiri et al.l 
(|2006bl ). Minor axis is found to be preferentially aligned with 
the direction to the centre of the void and major axis tends 
to be aligned with the orthogonal direction. These results 
are also robust with respect to the shape of the void. Us- 
ing an ellipsoidal void catalogue and measuring angles with 
respect to the radial and the normal direction we do not 
find any difference compared to the results obtained using 



spherical voids for the orientation of the main axes and an- 
gular momentum. Hence, we do not find that the shape of 
the void is meaningful for this purpose. We also analysed 
the strength of these alignments in several shells around the 
surface of the void. Comparing our results with those in 
iBrunino et al.l l|2007l ) we find a systematically weaker align- 
ment in the shape orientation of dark matter haloes. 

Regarding gas particles, we find that the angular mo- 
mentum of baryons is oriented with respect to the radial 
direction with a strength similar to dark matter, even con- 
sidering the inner parts of the halo. However, in the core of 
the halo shape alignment for gas particles is enhanced with 
respect to dark matter. Although this is a remarkable result 
for its observational implications, the effect of cooling and 
full scale star formation in these orientations remains to be 
investigated. 

In addition, we found a slight trend of the align- 
ment of angular moment um with the halo mass, as in 
lAragon-Calvo et al.l l|200it ). More massive haloes happen to 
be more aligned with the orthogonal direction. A possible 
explanation to this fact could be their recent formation, as 
pointed out by Aragon-Calvo et al.. (2007, ), which makes less 
massive haloes be more affected by non- linear interactions. 
It has also been suggested that an alternative explanation 
is that this effect might be due to th e preferred accretion 
along the filaments l|Shaw et al.ll2003 ). 

Evolution of this alignment in cosmological simulations 
have been studied by iHahn et al.l l|2007bf ). Unfortunately, 
it is very difficult to draw a conclusion about the evolution 
of the alignments from the present data. A deeper analysis 
of the redshift evolution of the strength of this signal is 
mandatory to elucidate its origin. 
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